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It is shown that the renormalisation constants of two quark operators can be accurately determined (to a 
precision of a few percent using 18 gluon configurations) using chiral Ward identities. A method for comput- 
ing renormalisation constants of generic composite operators without the use of lattice perturbation theory is 
proposed. 



1. Introduction 

In this talk I will report on attempts to de- 
velop techniques for the evaluation of renormali- 
sation constants of composite operators, without 
the use of lattice perturbation theory. These con- 
stants are needed to obtain physical results for 
quantities such as mesonic decay constants and 
weak and electromagnetic form-factors from lat- 
tice computations of the corresponding matrix el- 
ements. Although much of the discussion will be 
quite general, the specific applications considered 
below will be to matrix elements of the improved 
operators [1] 
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computed using the Sheihkoleslami and Wohlert 
(SW) "improved" fermion action [2]. The discreti- 
sation errors in these matrix elements are of 
O(asa) [1] (as compared to 0(a) errors in ma- 
trix elements computed using the standard Wil- 
son fermion action). 

In the next section I will review the use of chi- 
ral Ward identities for the determination of renor- 
malisation constants. [3]. The techniques used in 
sec. 2 are generalisations of those proposed in refs. 
[4] and [5] for simulations with the Wilson fermion 
action. In section 3 I will briefly describe the sta- 
tus of our attempts to compute the renormalisa- 
tion constants of general lattice operators. The 



numerical results presented in this talk were ob- 
tained from a simulation in which 18 independent 
gluon configurations were generated on a 16^ x 32 
lattice at /3 = 6.0. 

2. Chiral Ward Identities 

2.1. The Conserved Vector Current 

In the theory defined by the SW fermion ac- 
tion action there is a vector current, V^^ , which 
is both conserved and improved [6]. V^^ is con- 
served, and so has no perturbative corrections in 
forward matrix elements, i.e. Zyci = 1, and it 
is improved so that it has no 0(a) corrections 
even in non-forward matrix elements. Defining 
the "improved" vector current, , as that given 



in eq.(l) with F 
slant is given by 
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j^, its renormalisation con- 
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ZyV^ is the correctly normalised vector current 
[5,7]. In the simulation at /3 = 6.0, with \i) = \p) 
and I/) = |0), we find Zv = 0.824(2) at k = 
0.1425 (for which the pion has a mass of about 
900 MeV) and Zy = 0.833(1) at k = 0.1410 
(for which the pion has a mass of about 1100 
MeV). The statistical errors in the ratios (2) tend 
to be very small. These results can be com- 
pared to that obtained from perturbation theory 



2 



Zy = \- O.lOff^ ~ 0.83 [8], where the numeri- 
cal value was estimated using a boosted coupling 
ff2 = 1.7(6//3). 

2.2. The Renormalisation Constant of the 
Axial- Vector Current 

We now determine Zy and Z a by requiring 
that the currents Zy and ZaA^ satisfy the 
continuum chiral Ward identities. These iden- 
tities, derived from the axial flavour transforma- 
tions 6aiP(x) = ia'^j^(x)^j5ip(x) and SAipix) = 
ia{^(x)-ip(x) imply [3] 

2p J d'x J d^y{P''i{x)Ai'^{y)Vl'\Q)) = 

-i - pra^ ff^' I d'y{Ai'ny)AYm (3) 

and are the pseudoscalar density and axial- 
vector current respectively (defined by (1) with 
r = 7® or 7''7®) and 2p is given by 

2 ^ B,Y.y< Ai-fiyJy)P^'-fiQ,Q) > 

where ty ^ 0, so that the operators are sepa- 
rated (in practice, to avoid the effects of contact 
terms, we take ty to be greater than a few lat- 
tice spacings). A-'^ and f^^^ are the generators 
and structure constants of the flavour symme- 
try group respectively. The terms proportional 
to pra in eq.(3) come from contact terms arising 
from the use of the equations of motion in defining 
the improved operators in (1). By requiring that 
eq.(3) holds at each ty the combinations Zy IZ'a 
and XjZy (and hence Zy and Za) can be deter- 
mined. Zy and Za have no discretisation errors 
of 0{a). The result for Zy must, of course, be 
consistent with that obtained from the determi- 
nation using eq.(2). 

Since the Ward identity in eq.(3) involves con- 
tact terms, it may be thought that large momenta 
are present which might spoil improvement. The 
integration over x and y on the left hand side of 
eq.(3) ensures that this is not the case (at least 
in perturbation theory) [3]. 



Our best estimate for Za at /3 = 6.0 and 
K = 0.1425 is 1.09(3). Our preliminary value at 
K = 0.1410 is consistent with this (Za = 1.08(2)). 
The result from perturbation theory is Za = 
1 - 0.02ff2 ~ 0.97. The values of Zy obtained 
using three-point functions (with i> = p = 4) 
are consistent with those obtained from the two- 
point functions in subsec. 2.1 but have larger er- 
rors (Zy=0.85(3) at K=0.1425). 

The Ward Identity for the pseudoscalar and 
scalar densities gives the ratio of the correspond- 
ing renormalisation constants for which we find 
(at K = 0.1425) Zp/Zs = 0.60(2) (the result from 
perturbation theory is (l-0.2Sr/)/(l-0mr/) ~ 
0.68). 

2.3. Za Obtained Using Quark Green 
Functions 

We have also evaluated Za by studying the fol- 
lowing Ward identity for quark Green Functions: 

2/9 j d'^x j d^y{ua{y) {u{x)j5d{x)) d/siO)) = 

(^^-pra^ J d''y{^rA.y)m + u{y)u{Q)r')ai5 

where u and d are the up and down quark fields 
respectively. We have evaluated Za using eq.(5) 
in the Landau gauge and found Za = 1.14(8) at 
K = 0.1425 (Za = 1.08(6) at k = 0.1410). These 
results are in excellent agreement with those ob- 
tained using the Ward identity (3). It is pleasing 
to obtain results using quark Green's functions 
with reasonably small errors, and the role of Gri- 
bov copies in such calculations is being investi- 
gated. 

3. Towards Non-Perturbative Renormah- 
sation of General Composite Operators 

The technology of higher order perturbative 
calculations in continuum QGD is well advanced, 
two and three loop calculations are now routinely 
performed. In lattice perturbation theory the 
calculations are considerably more complicated, 
and the one-loop coefficients are frequently large. 
Lepage and Mackenzie have provided consider- 
able insights into the reasons for this [9] , and have 
suggested a procedure for an efficient reorgani- 
sation of the perturbation series, based on mean 
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field theory and the partial summation of tadpole 
graphs. Nevertheless, the ignorance of higher or- 
der perturbative coefficients in lattice perturba- 
tion theory for renormalisation constants of com- 
posite operators represents an important source 
of uncertainty in lattice computations. I conclude 
this talk by briefiy explaining an exploratory pro- 
gram we are pursuing in an attempt to avoid the 
need for lattice perturbation theory for matrix el- 
ements of general operators[10]. 

For the purposes of illustration I consider a 
multiplicatively renormalisable operator O. We 
propose to define the renormalised operator Or 
by fixing a renormalisation condition which can 
be imposed non-perturbatively on the lattice. For 
example it may be convenient to require that an 
S-matrix element of Or, taken between quark 
states at some fixed momentum, takes a specific 
value. In our first attempts to implement these 
ideas we evaluate ratios of the form 
S-\p)Go,{p)S-\p) 

s-Hp)GvMs-Kp) ^ ^ 

where S{P) is the quark propagator, 

Gob(p) = j d'^xd'^yexp{ip ■ x) {il){x)OB{y)i'{Q)) 

Ob is the bare operator (with the lattice cut-off) 
and Vr is the time component of the bare vector 
current. The normalisation by the Green function 
with the vector current is a convenient method of 
implementing the necessary wave function renor- 
malisation (i.e. multiplying by \JZ^ for each ex- 
ternal line). Our calculations are performed in 
the Landau gauge. The renormalised operator 
Or is defined by imposing the condition that the 
ratio (5), with Gqr in the numerator, takes a spe- 
cific value at momentum p, (in practice we trace 
the numerators and denominators of (5) with the 
corresponding 7 matrix). To obtain the matrix ele- 
ments in some other standard continuum scheme, 
e.g. the MS scheme, then requires a calculation 
in continuum perturbation theory only. Our ini- 
tial results have encouragingly small errors, here 
I only present a few sample preliminary results 
for finite quantities, from our simulation on the 
16^ X 32 lattice at /3 = 6.0, with p = (tt/S, 0) and 
K = 0.1425. For the ratios of (5) for the pseu- 
doscalar and scalar densities we find 3.3(4), for 



the axial current 0.62(5), and for Z^S~^Gvb^~^ 
we find 1.26(7), where for we take the coef- 
ficient of P47'* in the quark propagator. We are 
currently investigating the corresponding results 
in both lattice and continuum perturbation the- 
ory to gain some insights into our results, and 
their behaviour with momenta in particular [10]. 

An important question in our program is 
whether it is possible to make the momenta suf- 
ficiently small so that lattice artefacts are negli- 
gible {pa <ti 1), and yet sufficiently large that the 
perturbative calculations are valid (p/Aqcd ^ 
1). This is under investigation. 

Although the discussion presented here has 
been for quark bilinear operators, it is also likely 
to be very useful for four-quark operators relevant 
for weak decays, and for operators in the heavy 
quark effective theory. 
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